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Abstract 

Hydrogen  production  via  steam  electrolysis  may  involve  less  electrical  energy  consumption  than  conventional  low  temperature  water  electrolysis, 
reflecting  the  favourable  thermodynamics  and  kinetics  at  elevated  temperatures.  The  present  paper  reports  on  the  development  of  a  one-dimensional 
dynamic  model  of  a  cathode-supported  planar  intermediate  temperature  solid  oxide  electrolysis  cell  (SOEC)  stack  with  air  flow  introduced  through 
the  cells.  The  model,  which  consists  of  an  electrochemical  model,  two  mass  balances,  and  four  energy  balances,  is  here  employed  to  study  the 
prospect  of  the  stack  temperature  control  through  the  variation  of  the  air  flow  rate.  The  simulations  found  that  the  increase  in  the  air  flow  rate 
provides  enhanced  cooling  and  heating  during  exothermic  and  endothermic  operations,  respectively.  The  stack  behaviour  has  suggested  that  such  a 
convective  heat  transfer  between  the  cell  components  and  air  flow  would  allow  the  control  of  stack  temperature.  However,  only  a  small  dependence 
of  the  temperature  on  the  air  flow  rate  was  observed  for  a  stack  driven  at  conditions  near  thermoneutral  operation,  indicating  that  this  operating 
mode  should  be  avoided  from  a  control  perspective. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Hydrogen  is  regarded  as  a  leading  candidate  for  alternative 
future  fuels.  It  has  the  potential  to  address  the  environmental  and 
energy  security  issues  associated  with  fossil-derived  hydrocar¬ 
bon  fuels.  Among  many  hydrogen  production  methods,  water 
electrolysis  is  a  well-established  technique,  which  is  capable 
of  producing  carbon-free  hydrogen  if  used  in  conjunction  with 
renewable  or  nuclear  energy.  However,  water  electrolysis  has 
not  had  a  significant  commercial  penetration,  mainly  due  to  its 
high  electricity  consumption  and  associated  high  operating  cost 
[1]. 

Steam  electrolysis  at  elevated  temperatures  might  offer  a 
solution  by  consuming  less  electricity  than  is  required  at  ambi¬ 
ent  conditions  through  favourable  thermodynamics  and  kinetics 
[2],  Such  a  method  of  hydrogen  production  is  performed  using  a 
solid  oxide  electrolysis  cell  (SOEC),  which  can  be  seen  in  simple 
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terms  as  the  reverse  operation  of  a  solid  oxide  fuel  cell  (SOFC), 
allowing  the  opportunity  to  apply  recent  developments  in  SOFCs 
to  the  field.  An  SOEC  consists  of  a  three-layer  solid  structure 
(composed  of  porous  cathode,  electrolyte  and  porous  anode)  and 
an  interconnect  plate.  Steam  is  introduced  at  the  cathode  side  of 
the  solid  structure  where  it  is  reduced  to  hydrogen,  releasing 
oxide  ions  in  the  process.  The  oxide  ions  then  migrate  through 
the  electrolyte  to  the  anode  where  they  combine  to  form  oxy¬ 
gen  molecules,  releasing  electrons.  Although  the  use  of  proton 
conductors  has  also  been  investigated  [3-5],  yttria-stabilised  zir- 
conia  (YSZ),  which  is  an  oxide  ion  conductor,  is  generally  used 
for  the  electrolyte  in  SOECs.  Typical  materials  for  the  cathode 
are  nickel-YSZ  cermets  and  those  for  the  anode  are  perovskite 
oxides  such  as  lanthanum  manganite. 

In  SOFC  technologies,  there  is  increased  interest  in  interme¬ 
diate  temperature  SOFCs  (IT-SOFCs),  which  typically  operate 
between  823  and  1073  K,  to  allow  for  a  wider  range  of  materials, 
more  cost  effective  SOFC  fabrication  methods  [6]  and  increased 
durability.  The  reduction  in  operating  temperature  has  also  been 
seen  in  the  field  of  SOECs  in  the  limited  number  of  experi¬ 
mental  projects  reported  over  recent  decades  [7-9],  where  the 
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Nomenclature 

C,  concentration  of  gas  species  i  in  the  anode  gas 
streams  (molm-3) 

Co2 ,  Cn2  concentration  of  oxygen  and  nitrogen  in  the 
anode  gas  streams  (molm-3) 

C®  concentration  of  gas  species  i  at  the  stack  inlet 
(mol  m-3) 

C™  concentration  of  oxygen  at  the  triple  phase  bound¬ 

ary  (molm-3) 

f?eff,  anode  average  effective  diffusivity  coefficient  of  the 
anode  (m2  s-1) 

fl  inlet  flow  rate  of  the  anode  streams  (mol  s-1) 

F  Faraday’s  constant  (C  mol- 1 ) 

h\  height  of  the  anode  gas  channels  (m) 

j  local  current  density  (A  m-2) 

j  average  current  density  (A  m-2) 

jo,  anode  exchange  current  density  of  the  anode  (A  m-2) 

L  cell  length  (m) 

R  reaction  rate  (mol  s- 1  m-2) 

ill  gas  constant  (J  mol-1  K-1) 

t  time  (s) 

7a,  7’s  temperature  of  the  anode  gas  streams  and  solid 
structure  (K) 

ua  velocity  of  the  anode  gas  streams  (ms-1) 

W  cell  width  (m) 

x  axial  coordinate  (m) 

y™le  mole  fraction  of  oxygen  in  the  anode  gas  streams 
Greek  letters 

a  transfer  coefficient 

?7act,  anode  activation  overpotential  losses  at  the  anode  (V) 
r)conc  anode  concentration  overpotential  losses  at  the  anode 
(V) 

Vi  stoichiometric  coefficient  of  gas  species  i 

vq2  stoichiometric  coefficient  of  oxygen 

tanode  thickness  of  the  anode  (m) 
i//  air  ratio 


temperature  at  which  the  cells  are  tested  has  decreased  from 
1273  to  1073  K.  Initially  SOECs  were  developed  as  tubular  cell 
structures  to  avoid  sealing  problems,  which  is  a  major  issue  in 
segregating  the  produced  H2  and  O2  in  planar  cells,  particularly 
over  multiple  thermal  cycles  [10],  Today,  despite  such  sealing 
difficulties,  the  trend  in  cell  design  is  to  employ  planar  structures. 
A  planar  structure  permits  high  packing  density  and  significantly 
smaller  hot  volume  in  the  system  than  that  allowed  by  a  tubu¬ 
lar  design  [11],  together  with  lower  manufacturing  costs  and 
shorter  current  paths,  reducing  the  ohmic  cell  resistance  [12],  A 
planar  SOEC  can  be  categorised  as  being  either  an  electrolyte- 
supported  or  electrode-supported  design.  The  former  employs 
the  electrolyte  as  the  support  structure  and  is  suitable  for  high 
temperature  operation  in  which  the,  often  large,  ohmic  resis¬ 
tance  associated  with  a  thick  electrolyte  can  be  reduced.  In  an 
electrode-supported  cell,  on  the  other  hand,  one  of  the  electrodes 


an  inlet  temperature  of  1023  K. 


is  the  thickest  part  of  the  solid  structure,  supporting  a  thick  film 
electrolyte  typically  10  p,m  or  so  in  thickness.  Such  a  design 
has  been  developed  in  order  to  minimise  ohmic  resistances  in 
SOFCs  operating  at  intermediate  temperatures,  and  may  also  be 
applied  to  intermediate  temperature  SOECs  (IT-SOECs). 

To  ensure  a  sufficient  rate  of  H2  production,  an  SOEC  sys¬ 
tem  must  consist  of  several  repeating  cells  assembled  in  stacks. 
An  SOEC  stack  can  function  in  either  exothermic,  endother¬ 
mic  or  thermoneutral  operating  modes.  During  the  operation  of 
an  SOEC  system,  heat  is  generated  in  the  stack  as  the  result  of 
irreversible  losses  due  to  ohmic  resistance  and  electrode  overpo¬ 
tentials.  These  depend  on  operating  conditions  such  as  the  stack 
temperature  and  average  current  density,  as  well  as  materials 
selection  and  cell  and  stack  geometry.  Fig.  1  shows  examples  of 
stack  temperature  distribution,  produced  by  the  one-dimensional 
model  of  a  co-flow  IT-SOEC  stack  described  in  a  previous  pub¬ 
lication  [13].  The  modelled  stack  considers  the  cells  with  the 
flow  of  H2  and  H2O  on  the  cathode  side,  and  of  pure  O2  on  the 
anode  side  of  the  solid  structure.  As  the  mixture  of  H2  and  H2O 
travels  along  the  stack  (from  left  to  right  in  the  figure),  H2O  is 
consumed  by  the  reaction  while  FE  and  O2  are  produced.  At 
the  inlet  to  the  cathode  gas  channels,  10mol%  H2  is  assumed 
sufficient  in  preventing  the  oxidation  of  the  electrode  material 
[7],  while  the  inlet  to  the  anode  gas  channels  are  sealed  to  allow 
the  collection  of  pure  O2  from  the  outlet.  An  inlet  temperature 
of  1023  K,  an  operating  pressure  of  0. 1  MPa  and  a  steam  util¬ 
isation  factor  of  80%  have  been  selected.  For  the  steady  state 
simulation  with  an  average  current  density  of  7000  Am-2,  the 
heat  generated  via  irreversible  losses  exceeds  the  thermal  energy 
consumed  by  the  endothermic  electrolysis  reaction.  Such  a  con¬ 
dition  is  referred  to  as  exothermic  stack  operation  in  which  the 
temperature  increases  along  the  stack  due  to  the  heat  accumula¬ 
tion,  as  indicated  by  the  figure.  Although  exothermic  operation 
provides  the  opportunity  to  bring  the  inlet  gas  streams  to  the 
operating  temperature  entirely  through  the  recovery  of  the  heat 


356 


J.  Udagawa  el  al.  /Journal  of  Power  Sources  180  (2008)  354-364 


from  the  hotter  outlet  streams,  such  an  operating  mode  is  char¬ 
acterised  by  an  increased  electrical  energy  consumption  by  the 
stack.  Conversely,  for  the  simulation  with  an  average  current 
density  of  5000  A  m-2,  the  heat  generated  via  irreversible  losses 
is  predicted  to  be  smaller  than  the  thermal  energy  consumed 
by  the  reaction.  Such  a  situation  corresponds  to  endothermic 
stack  operation  in  which  the  temperature  decreases  as  the  reac¬ 
tion  proceeds  along  the  stack.  Although  endothermic  operation 
allows  the  stack  to  function  with  a  lower  electrical  energy  con¬ 
sumption,  the  rate  of  H2  production  per  unit  stack  area  is  less 
than  that  arising  from  exothermic  operation,  evident  from  the 
lower  average  current  density.  Furthermore,  during  endother¬ 
mic  operation  it  is  not  possible  to  bring  the  inlet  streams  to  the 
operating  temperature  entirely  by  heat  recovery  from  the  out¬ 
let  streams.  The  implication  is  that  an  external  heat  source  is 
required  to  supply  the  thermal  energy  necessary  to  elevate  the 
inlet  stream  temperature  to  the  required  level.  Therefore,  in  the 
design  of  SOEC  systems,  there  is  a  trade-off  between  the  ben¬ 
efits  of  partially  replacing  electrical  energy  by  less  expensive 
thermal  energy,  against  the  extra  costs  involved  in  engineering 
the  transfer  of  this  thermal  energy  from  external  sources.  Finally, 
thermoneutral  stack  operation  occurs  when  the  thermal  energy 
consumed  by  the  reaction  is  precisely  matched  by  the  heat  gen¬ 
erated  via  irreversible  losses.  During  this  mode  of  operation, 
although  the  outlet  streams  carry  the  same  amount  of  thermal 
energy  as  the  inlet  streams,  the  energy  losses  associated  with  the 
heat  recovery  process  mean  that  an  external  heat  source  and  cor¬ 
responding  equipment  costs  would  still  be  necessary,  as  in  the 
case  of  endothermic  operation.  The  electrical  energy  consump¬ 
tion  of  the  stack  during  thermoneutral  operation  is,  however, 
always  greater  than  that  during  endothermic  operation.  There¬ 
fore,  in  general,  thermoneutral  operation  is  not  recommended 
from  an  economic  perspective  [14], 

In  a  previous  publication,  a  one-dimensional  distributed 
dynamic  model  of  a  cathode-supported  planar  IT-SOEC  stack 
has  been  presented  [13].  The  steady  state  simulation  of  the  model 
predicted  an  electricity  consumption  significantly  less  than  that 
of  low  temperature  stacks  commercially  available  today.  How¬ 
ever,  the  dependence  of  the  stack  temperature  distribution  on  the 
average  current  density,  as  can  be  seen  in  Fig.  1,  calls  for  strict 
temperature  control  if  such  a  stack  is  to  be  successfully  built  and 
used,  especially  in  dynamic  operation.  Although  SOEC  systems 
are  often  considered  for  large  scale  steady  state  operation  in 
which  the  input  power  source  may  be  a  nuclear  reactor,  tem¬ 
perature  control  becomes  essential  if  the  systems  are  to  be  used 
in  dynamic  operation  with  intermittent  electrical  power  sources 
such  as  wind  turbines  or  photovoltaic  cells.  To  prevent  the  frac¬ 
ture  of  delicate  stack  components  during  dynamic  operation, 
significant  thermal  excursions  in  the  stack  need  to  be  avoided 
by  the  implementation  of  an  effective  control  strategy.  Such  a 
control  strategy  would  also  ensure  a  constant  operating  mode  of 
the  stack,  which  is  an  important  constraint  as  the  system  require¬ 
ments  vary  depending  on  whether  the  stack  is  employed  in  an 
exothermic,  endothermic  or  thermoneutral  mode.  One  possible 
solution  to  achieving  the  control  requirements  would  be  to  mod¬ 
ify  the  SOEC  design  such  that  one  more  degree  of  freedom  is 
available  to  control  the  stack  temperature.  As  can  be  seen  in  [15], 


the  temperature  of  SOFC  stacks  is  often  controlled  by  varying 
the  air  flow  through  the  cells.  Such  a  strategy  may  potentially 
be  applied  in  the  case  of  an  SOEC  stack  where  the  introduction 
of  air  to  the  anode  gas  channels  would,  in  addition  to  control¬ 
ling  the  stack  temperature,  dilute  the  O2  by-product.  Although 
this  would  remove  the  possibility  of  generating  extra  revenues 
through  the  sale  of  pure  O2,  it  would  also  reduce  the  risk  of  high 
temperature  corrosion,  which  might  otherwise  be  caused  by  the 
exposure  of  metallic  stack  and  balance  of  plant  components  to 
such  a  highly  oxidising  atmosphere  at  elevated  temperatures 
[16]. 

Note  that  the  introduction  of  the  air  flow  to  the  anode  channels 
would  involve  extra  investment  costs  for  the  compressors  and 
heat  exchangers  as  well  as  the  energy  costs  associated  with  the 
operation  of  such  additional  system  components.  These  costs 
would  be  a  function  of  the  operating  conditions,  system  scales 
as  well  as  the  origin  of  the  input  energy.  In  the  development  of 
the  control  strategy,  it  is  important  to  evaluate  such  extra  costs 
and  their  magnitude  with  respect  to  the  overall  cost  of  the  H2 
production  as  well  as  to  investigate  the  technical  feasibility  of 
the  control  strategy. 

A  mathematical  model  is  an  important  design  tool  for  devices 
such  as  an  SOEC,  which  are  still  in  the  development  stage.  It 
allows  the  prediction  of  the  behaviour  of  the  device  under  dif¬ 
ferent  process  conditions  and  assists  in  the  optimisation  of  its 
performance.  In  particular,  such  a  model  is  essential  in  under¬ 
standing  the  response  of  an  electrolyser  under  steady  electrical 
power  input,  such  as  those  from  the  grid  or  nuclear  energy, 
as  well  as  under  an  intermittent  renewable  electrical  power 
input.  However,  no  studies  on  the  control  issues  or  the  dynamic 
response  of  an  SOEC  can  be  found  in  the  literature.  The  present 
paper  reports  on  the  development  of  a  dynamic  model  of  a  con¬ 
ceptual  IT-SOEC  stack  with  the  air  flow  introduced  through 
the  cells,  produced  within  the  gPROMS  modelling  environ¬ 
ment.  This  paper  introduces  the  IT-SOEC  stack  model  and  the 
results  describing  the  simulated  steady  state  behaviour  of  such 
an  electrolyser.  The  prospect  of  controlling  temperature  tran¬ 
sitions  caused  by  the  time  dependent  demand  for  hydrogen 
from  the  IT-SOEC  is  discussed.  While  the  developed  model  has 
dynamic  capability,  only  the  steady  state  results  are  analysed 
here.  The  dynamic  response  will  be  presented  in  a  subsequent 
publication. 

2.  IT-SOEC  mathematical  model 

Although  the  number  of  modelling  activities  on  electroly¬ 
sis  at  elevated  temperatures  has  increased  over  recent  years 
[13,17-30],  the  mathematical  modelling  of  SOEC  stacks  is  still 
not  an  active  research  area.  Today,  only  limited  modelling  stud¬ 
ies  attempting  to  describe  SOEC  stack  behaviour  can  be  found 
[13,25-30]  and,  to  the  authors’  knowledge,  the  issue  of  tem¬ 
perature  control  within  an  SOEC  stack  has  not  been  the  focus 
of  any  published  work.  Here,  a  previously  developed  dynamic 
model  of  a  cathode-supported  planar  IT-SOEC  stack  [13]  has 
been  modified  to  introduce  air  flow  as  a  method  of  control¬ 
ling  stack  temperature.  Important  changes  implemented  to  the 
previously  developed  model  are  described  in  Sections  2. 1-2.4. 
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Fig.  2.  Schematic  view  of  a  planar  SOEC  stack  with  the  air  flow  through  the 
anode  channels. 


concentration  overpotential  was  included,  anode  concentration 
overpotential  was  not  taken  into  account,  on  the  basis  that  the 
difference  between  the  CH  concentration  at  the  triple  phase 
boundaries  TPBs  and  that  in  the  bulk  streams  was  assumed 
negligible.  However,  when  air  is  introduced  through  the  anode 
channels,  such  an  assumption  no  longer  necessarily  applies 
and  it  becomes  important  to  include  the  anode  concentration 
overpotential  in  the  model.  Eq.  (1)  describes  the  anode  concen¬ 
tration  overpotential  for  the  IT-SOEC  stack  with  air-fed  anode 
channels.  Note  that  the  temperature  of  the  gas  mixture  at  the 
TPBs  is  represented  by  the  solid  structure  temperature.  The  O2 
concentration  at  the  TPBs  is  determined  for  one-dimensional 
self-diffusion  of  O2  in  Eq.  (2).  Z)eff,  anode  represents  the  aver¬ 
age  effective  diffusivity  coefficient  in  the  anode,  which  is  taken 
as  13.7  x  10-6  m-2  s_1 ,  considering  a  binary  gas  mixture  of  O2 
and  N2  [6] .  A  method  for  expressing  the  mass  diffusion  activities 
in  an  SOFC  [32,33]  has  been  adopted  for  an  SOEC  in  deriving 
Eq.  (2). 


In  an  SOEC  system,  several  repeating  cells  are  assembled  in 
stacks  to  support  a  sufficient  rate  of  H2  production.  However,  the 
models  of  such  stacks  are  usually  constructed  for  the  smallest 
unit  cell,  which  is  assumed  to  describe  the  response  of  the  whole 
stack,  subject  to  the  use  of  adequate  boundary  conditions.  Here, 
the  modelled  unit  cell  is  considered  to  be  in  the  centre  of  a  large 
stack  such  that  end  effects  can  be  neglected.  Although  intercon¬ 
nects  normally  provide  the  gas  flow  channels,  above  and  below 
the  solid  structure,  the  effect  of  individual  passages  is  neglected. 
The  pressure  drop  along  the  gas  channels  is  also  assumed  to  be 
negligible  at  the  operating  pressure  of  0.1  MPa,  and  for  the  flow 
rates  considered.  For  the  purpose  of  developing  the  model,  the 
unit  cell  is  considered  to  be  composed  of  four  components,  the 
cathode  and  anode  gas  streams,  solid  structure  and  intercon¬ 
nect.  Fig.  2  shows  the  schematic  view  of  such  a  unit  cell  with 
the  air  flow  through  the  anode  gas  channel.  The  model  consists 
of  an  electrochemical  model,  mass  balances  for  the  cathode  and 
anode  streams,  and  energy  balances  for  the  cathode  and  anode 
streams,  solid  structure  and  interconnect.  The  properties  of  the 
gas  streams,  solid  structure  and  interconnect  are  assumed  to  be 
constant  in  order  to  achieve  a  more  computationally  tractable 
model.  This  simplification  has  been  assessed  in  equivalent  stud¬ 
ies  of  IT-SOFC  stacks  in  which  the  stack  model  employing  the 
constant  flow  properties,  determined  at  the  stack  inlet,  and  that 
employing  the  variable  flow  properties  have  been  compared  [31]. 
The  study  showed  that  the  simplification  is  justified  provided 
the  average  current  density  is  not  high.  For  the  IT-SOEC  stack 
model  here,  to  minimise  the  errors  involved  in  using  the  constant 
properties  even  at  relatively  high  average  current  densities,  the 
gas  properties  are  averaged  over  the  entire  stack  length,  rather 
than  using  those  determined  at  the  stack  inlet.  Finally,  ideal  gas 
behaviour  is  assumed  for  the  cathode  and  anode  streams. 

2.1.  Anode  concentration  overpotential  losses 

Concentration  overpotentials  occur  as  mass  transport  of  the 
gas  species  involved  in  the  reaction  becomes  increasingly  rate 
limiting.  In  the  previously  developed  model,  whilst  cathode 
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2.2.  Anode  activation  overpotential  losses 

Activation  overpotentials  are  related  to  the  chemical  kinetics 
of  the  reactions.  They  are  classically  determined  through  the 
Butler- Volmer  equation.  In  the  previously  developed  model,  as 
the  difference  between  the  O2  concentration  at  the  TPBs  and 
that  in  the  bulk  stream  was  assumed  negligible,  the  simplified 
form  of  the  Butler- Volmer  equation  was  applied  in  predicting 
the  anode  activation  overpotential.  Here,  the  simplified  form 
has  been  replaced  by  the  extended  form  (3)  to  account  for  the 
difference  in  the  O2  concentration  created  between  the  TPBs 
and  bulk  stream  by  the  introduction  of  the  air  flow  through  the 
anode  channels. 

j(x)  =  JO.anodeO)  |exp  anode 

CofW  [-2  aF  1 

Co20)  L'K7sW  J 

2.3.  Anode  stream  mass  balance 

Although  a  mass  balance  can  track  the  changes  in  the  stream 
compositions,  this  was  not  applied  to  the  anode  streams  in  the 
previously  developed  model,  as  the  streams  of  pure  O2  would 
involve  no  such  changes  along  the  stack.  Therefore,  for  the  IT- 
SOEC  stack  with  air-fed  anode  channels,  an  anode  stream  mass 
balance  must  be  introduced  in  the  model  because  the  compo¬ 
sitions  of  the  anode  streams  evolve  as  the  mixture  of  O2  and 
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N2  travels  towards  the  outlet.  Here,  the  anode  stream  mass  bal¬ 
ance  (4)  and  the  boundary  conditions  at  the  inlet  (5)  predict  the 
O2  and  N2  concentrations  at  each  location  along  the  stack.  The 
anode  stream  velocity  is  assumed  constant  and  determined  from 
the  outlet  flow  rate  of  the  streams. 

ViR(x),  i  e  {02,  N2}  (4) 

3 1  ox  h\ 

C1(0)  =  Cf,  i  e  {02,  N2}  (5) 


2.4.  Air  ratio 


The  air  ratio  reflects  the  inlet  flow  rate  of  air  in  relation  to 
the  rate  of  reaction.  It  is  here  defined  as  the  ratio  between  the 
moles  of  O2  contained  in  the  inlet  air  flow  to  that  produced 
in  the  unit  cell,  per  unit  time.  In  spite  of  the  slightly  different 
physical  definition,  the  mathematical  description  of  the  air  ratio, 
provided  in  Eq.  (6)  for  an  SOEC,  is  equivalent  to  that  of  an  SOFC 
[6].  As  will  be  discussed  in  the  following  sections,  it  would  be 
possible  to  provide  the  temperature  control  for  an  SOEC  stack 
through  the  manipulation  of  the  air  ratio.  Assuming  a  minimum 
of  50  mol%  N2  in  the  anode  streams  at  the  stack  outlet  is  required 
to  limit  the  corrosion  of  metallic  components,  the  lower  bound 
for  the  air  ratio  can  be  selected  to  be  as  small  as  0.4.  The  upper 
bound  is  assumed  to  be  14,  namely  the  maximum  air  flow  rate 
which  can  be  supplied  without  incurring  significant  additional 
energy  costs  [15].  Note  that  for  a  high  air  ratio,  the  anode  stream 
velocity  is  dominated  by  the  large  flow  rate  of  the  inlet  air.  Under 
such  conditions,  the  increase  in  the  anode  stream  mass  along  the 
stack,  which  results  from  the  accumulation  of  the  product  O2,  has 
negligible  influence  over  the  stream  velocity.  However  as  the  air 
ratio  is  lowered  and  the  flow  rate  of  the  inlet  air  is  consequently 
reduced,  such  a  rising  stream  mass  becomes  increasingly  more 
important  in  accelerating  the  stream  towards  the  outlet.  It  is, 
therefore,  important  to  point  out  that  any  inaccuracies  incurred  in 
assuming  a  uniform  anode  stream  velocity  in  the  present  model 
can  become  progressively  larger  as  the  air  ratio  is  reduced. 


o  jp-fO  ,,mole 

2/7Ay02 

]LWvo2 
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3.  Simulation  results  and  discussion 


With  the  inclusion  of  Eqs.  (1),  (2)  and  (4)-(6)  into  the  pre¬ 
viously  developed  model  and  the  replacement  of  the  simplified 
Butler-Volmer  equation  for  the  anode  activation  overpotential 
by  the  extended  form  (3),  the  system  of  partial  differential 
and  algebraic  equations  for  the  model  of  an  IT-SOEC  stack 
with  air-fed  anode  channels  is  solved  via  the  finite  difference 
method  using  gPROMS  Model  Builder  3.0.3  [34],  The  cath¬ 
ode  and  anode  inlet  compositions  are  assumed  to  be  10mol% 
H2/90mol%  H20  and  21mol%  O2/79  mol%  N2,  respectively. 
The  steam  utilisation  factor  is  selected  to  be  80%.  The  stack 
geometry  and  material  properties  are  found  in  the  previous  pub¬ 
lication.  Section  3.1  introduces  the  steady  state  performance 
of  the  IT-SOEC  stack  with  air-fed  anode  channels,  comparing 
the  reversible  potential  and  irreversible  losses  to  those  of  the 


Dimensionless  axial  position 

Fig.  3.  Reversible  potential  and  irreversible  losses  along  the  stack  for  an  average 
current  density  of  7000  A  m-2 ,  an  inlet  temperature  of  1 023  K  and  an  air  ratio 
of  7. 

IT-SOEC  stack  with  pure  O2  anode  streams.  The  proposed  tem¬ 
perature  control  strategy  through  the  manipulation  of  the  air  ratio 
is  explored  in  Section  3.2  in  which  the  stack  temperature  distri¬ 
butions  are  presented  for  different  values  of  the  air  ratio.  Section 
3.3  discusses  the  links  between  such  temperature  distributions 
and  the  local  current  density  along  the  stack. 

3.1.  Steady  state  performance  of  IT-SOEC  stack  with 
air-fed  anode  channels 

Figs.  3  and  4  present  the  contributions  of  the  reversible  poten¬ 
tial  and  irreversible  losses  on  the  cell  potential  for  the  IT-SOEC 
stack  with  air-fed  anode  channels  at  average  current  densities  of 
7000  and  5000  Am-2,  representing  exothermic  and  endother- 


current  density  of  5000 Am  ,  an  inlet  temperature  of  1023 K  and  an  air  ratio 
of  7. 
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mic  operation,  respectively.  In  both  cases,  the  inlet  temperature 
of  the  anode  streams  is  set  equivalent  to  that  of  the  cathode 
streams  at  1023  K  to  prevent  significant  changes  in  the  tempera¬ 
ture  along  the  stack  near  the  inlet,  and  the  air  ratio  has  been  fixed 
at  7.  Note  that,  as  observed  in  Fig.  1,  the  inlet  temperature  cor¬ 
responds  to  the  lowest  temperature  along  the  exothermic  stack 
and  highest  temperature  along  the  endothermic  stack.  The  intro¬ 
duction  of  the  air  flow  at  such  a  temperature,  therefore,  results 
in  the  convective  heat  transfer  between  the  cell  components  and 
air  flow,  providing  cooling  and  heating  to  the  exothermic  and 
endothermic  stacks,  respectively.  Furthermore,  the  introduction 
of  the  air  flow  causes  a  decrease  in  the  O2  partial  pressure  in 
the  anode  channels,  creating  the  difference  between  the  O2  con¬ 
centration  at  the  TPBs  and  that  in  the  bulk  streams.  Although 
this  is  accounted  for  in  the  prediction  of  the  anode  concentration 
overpotential  by  Eq.  (1),  both  Figs.  3  and  4  indicate  that  such  an 
irreversible  loss  is  negligible  for  the  conditions  employed.  As 
was  the  case  in  the  IT-SOEC  stack  with  pure  O2  anode  streams 
[13],  the  activation  overpotentials  are  the  predominant  source  of 
the  irreversible  losses. 

In  Fig.  3,  the  reversible  potential  for  the  exothermic  operation 
of  the  IT-SOEC  stack  with  air-fed  anode  channels  is  predicted 
to  be  around  0.027  V  lower  than  that  reported  in  the  previous 
publication  for  the  IT-SOEC  stack  with  pure  O2  anode  streams 
under  the  same  conditions.  As  illustrated  by  the  Nemst  equation, 
such  a  reduction  in  the  reversible  potential  is  a  consequence  of 
the  decrease  in  the  O2  partial  pressure  in  the  anode  channels. 
Although  the  reduced  stack  temperature,  which  results  from  the 
convective  cooling  supplied  by  the  air  flow,  tends  to  increase 
the  reversible  potential,  such  an  effect  is  overwhelmed  by  the 
influence  of  the  O2  partial  pressure  change  here.  The  steady 
state  simulations  suggest  that  the  introduction  of  the  air  flow,  at 
the  air  ratio  of  7,  reduces  the  average  stack  temperature  from 
1038  to  1027  K  and  the  average  O2  partial  pressure  from  0. 10  to 
0.02  MPa.  By  only  accounting  for  the  stack  temperature  change 
in  the  Nemst  equation,  such  a  temperature  reduction  alone  is 
expected  to  increase  the  average  reversible  potential  by  around 
0.005  V.  On  the  other  hand,  the  decrease  in  the  average  reversible 
potential  which  would  result  only  from  the  effect  of  the  O2  par¬ 
tial  pressure  reduction  is  estimated  to  be  as  much  as  0.034  V.  The 
impact  of  the  reduced  O2  partial  pressure  is  therefore  approx¬ 
imated  to  be  7  times  that  of  the  reduced  stack  temperature  for 
the  conditions  employed.  However,  in  spite  of  the  decreased 
reversible  potential,  the  cell  potential  shown  in  Fig.  3  is  almost 
equivalent  to  that  of  the  IT-SOEC  stack  with  pure  O2  anode 
streams.  This  is  due  to  the  compensating  increase  in  the  ohmic 
losses  and  activation  overpotentials,  of  around  0.034  V  on  aver¬ 
age  along  the  stack,  caused  by  the  fall  in  the  stack  temperature. 

Similarly  to  Fig.  3,  Fig.  4  shows  the  reversible  potential  for  the 
IT-SOEC  stack  with  air-fed  anode  channels  to  be  lower  than  that 
reported  in  the  previous  publication  for  the  IT-SOEC  stack  with 
pure  O2  anode  streams  during  endothermic  operation.  However, 
noticeable  differences  in  the  irreversible  losses  are  not  observed 
here  because  the  temperature  of  the  endothermic  stack  has  not 
been  altered  significantly  by  the  introduction  of  the  air  flow.  As 
will  be  explained  in  Section  3.2,  the  reduction  in  the  O2  partial 
pressure,  associated  with  the  dilution  of  the  anode  streams  by  the 


air  flow,  increases  the  thermal  energy  consumed  by  the  electrol¬ 
ysis  reaction  per  mole  of  product  H2.  Such  an  increased  thermal 
energy  consumption  results  in  the  increased  amount  of  heat 
removed  from  the  stack  and  has  an  effect  of  decreasing  the  stack 
temperature.  Although  this  effect  complements  the  convective 
cooling  supplied  by  the  air  flow  during  exothermic  operation, 
it  competes  against  the  convective  heating  during  endothermic 
operation.  Therefore,  the  temperature  of  the  endothermic  stack 
is  not  altered  to  the  same  extent  as  that  of  the  exothermic  stack 
by  the  introduction  of  the  air  flow. 

3.2.  Stack  temperature  control  through  air  ratio 
manipulation 

The  temperature  distribution  along  the  IT-SOEC  stack  with 
air-fed  anode  channels  operated  at  average  current  densities  of 
7000  and  5000  A  m-2  are  respectively  illustrated  in  Figs.  5  and  6. 
The  air  ratio  is  varied  between  0.4  and  14  while  the  inlet  tem¬ 
perature  is  maintained  at  1023  K.  As  can  be  seen  from  both 
figures,  the  increase  in  the  air  ratio  causes  the  stack  temperature 
to  become  more  uniformly  distributed,  approaching  the  temper¬ 
ature  of  the  inlet  streams.  Such  a  stack  behaviour  suggests  that 
the  increase  in  the  air  ratio  provides  enhanced  cooling  for  the 
stack  during  exothermic  operation,  and  enhanced  heating  during 
endothermic  operation,  via  the  increase  in  the  convective  heat 
transfer  between  the  cell  components  and  air  flow.  As  detailed 
previously  [13],  an  increase  in  average  current  density  causes 
the  irreversible  losses  to  increase.  During  exothermic  operation, 
the  consequent  rise  in  stack  temperature  may  be  counteracted 
by  increasing  the  air  ratio.  Conversely,  if  the  average  current 
density  is  reduced  during  exothermic  operation,  the  fall  in  stack 
temperature  may  be  avoided  by  decreasing  the  air  ratio.  Dur¬ 
ing  endothermic  operation,  on  the  other  hand,  the  rise  in  the 
stack  temperature  associated  with  the  increased  average  current 
density  may  be  balanced  by  a  decrease  in  the  air  ratio,  and  the 
fall  in  the  stack  temperature  corresponding  to  a  reduction  in  the 


Fig.  5 .  Cathode  stream  temperature  along  the  stack  for  an  average  current  density 
of  7000  A  m-2,  an  inlet  temperature  of  1023  K  and  air  ratios  of  0.4,  7  and  14. 
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Dimensionless  axial  position 

Fig.  6.  Cathode  stream  temperature  along  the  stack  for  an  average  current  density 
of  5000  A  m-2,  an  inlet  temperature  of  1023  K  and  air  ratios  of  0.4,  7  and  14. 

average  current  density  may  be  compensated  by  an  increased  air 
ratio. 

In  spite  of  such  a  concept  for  the  temperature  control,  closer 
observation  of  Figs.  5  and  6  reveals  that,  even  at  the  outlet,  the 
stack  temperature  only  varies  over  the  entire  range  of  air  ratios 
studied  by  around  13  K  for  the  exothermic  operation,  and  less 
than  7  K  for  the  endothermic  operation.  Here,  the  weak  depen¬ 
dence  of  the  stack  temperature  on  the  air  ratio  is  explained  by  the 
fact  that  the  sets  of  operating  conditions  used  in  producing  these 
results  are  close  to  thermoneutral  operation.  As  illustrated  later 
in  this  section,  if  a  stack  is  instead  operated  at  highly  exother¬ 
mic  or  endothermic  conditions,  the  magnitude  of  the  temperature 
gradient  would  be  large  along  the  stack  and  considerable  tem¬ 
perature  difference  would  exist  between  the  inlet  air  and  rest  of 
the  stack.  Under  such  conditions,  a  significant  transfer  of  heat 
between  the  cell  components  and  air  flow  would  occur,  and  the 
regulation  of  such  a  heat  transfer,  by  the  changes  in  the  air  flow 
rate,  should  result  in  a  satisfactory  temperature  control.  On  the 
other  hand,  if  a  stack  is  operated  close  to  the  thermoneutral  point, 
the  magnitude  of  the  temperature  gradient  would  be  small  and 
the  temperature  of  the  inlet  air  flow  would  be  similar  to  that  of 
the  rest  of  the  stack.  As  a  consequence,  only  a  modest  transfer 
of  heat  between  the  cell  components  and  air  flow  would  occur. 
The  regulation  of  such  a  minor  heat  transfer,  by  the  changes  in 
the  air  flow  rate,  is  not  expected  to  create  a  significant  influence 
over  the  stack  temperature.  Therefore  the  air  ratio  manipulation 
is  not  considered  to  provide  a  satisfactory  temperature  control 
for  the  conditions  presented  in  Figs.  5  and  6. 

Furthermore,  if  a  stack  is  operated  close  to  the  thermoneutral 
point,  there  is  an  increased  risk  of  encountering  the  tran¬ 
sition  between  exothermic  and  endothermic  operation  even 
with  the  aid  of  the  control  strategy.  For  the  conditions  pre¬ 
sented  in  Figs.  5  and  6,  the  decrease  in  the  average  current 
density  from  7000  to  5000 Am-2  would  result  in  the  stack 
to  switch  from  exothermic  to  endothermic  mode.  The  selec¬ 
tion  of  the  operating  conditions  such  that  the  stack  is  driven 


in  more  exothermic  or  endothermic  conditions  would  prevent 
such  transitions,  avoiding  the  necessary  changes  in  the  system 
configurations. 

In  addition  to  the  convective  transfer  of  heat,  the  changes  in 
the  air  flow  has  another  effect  on  the  stack  temperature  through 
the  changes  in  O2  partial  pressure.  The  thermodynamics  indi¬ 
cate  that,  assuming  ideal  gas  behaviour  of  H2,  H2O  and  O2,  the 
enthalpy  change  involved  during  the  electrolysis  reaction  is  inde¬ 
pendent  of  the  partial  pressure  of  the  gases  [35].  This  implies 
that  the  total  energy  consumed  by  the  reaction,  per  mole  of  prod¬ 
uct  H2,  is  not  altered  by  the  fall  in  the  O2  partial  pressure  caused 
by  the  dilution  of  the  anode  streams  by  the  air  flow.  In  con¬ 
trast,  the  electrical  energy  consumed  by  the  reaction,  per  mole 
of  product  H2,  is  decreased  by  the  fall  in  the  O2  partial  pressure. 
This  is  evident  in  the  previous  sections  where  the  introduction 
of  the  air  flow  to  the  anode  channels  caused  the  reduction  in  the 
reversible  potential  which  is  directly  proportional  to  the  electri¬ 
cal  energy  consumed  by  the  reaction.  During  electrolysis,  such 
a  decrease  in  the  electrical  energy  consumption  is  compensated 
by  an  increase  in  the  thermal  energy  consumption,  estimated 
in  Fig.  7  as  a  function  of  air  ratio.  The  figure  has  been  pro¬ 
duced  using  only  the  electrochemical  model,  without  mass  and 
energy  balances,  predicting  the  stack  behaviour  at  one  location 
along  the  stack,  which  is  assumed  to  be  at  1023  K  here.  The 
cathode  stream  composition,  averaged  over  the  stack  length,  is 
taken  as  46  mol%  H2/54  mol%  H2O.  The  average  anode  stream 
compositions  are  determined  for  each  air  ratio  values.  In  the  fig¬ 
ure,  the  thermal  energy  consumed  increases  from  around  66.9 
to  68.6  kJ  mol-1  H2  between  the  air  ratios  of  0.4  and  4  while 
a  significantly  smaller  rise  is  observed  above  the  air  ratio  of  4. 
This  is  because  the  reduction  in  the  O2  partial  pressure  involved 
when  the  air  ratio  is  increased,  for  example,  from  1  to  2  is  consid¬ 
erably  larger  than  that  involved  when  the  air  ratio  is  increased 
from  7  to  8.  During  exothermic  operation,  the  increase  in  the 
thermal  energy  consumption,  associated  with  the  rise  in  the  air 
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Fig.  7.  Thermal  energy  consumed  by  the  reaction,  per  mole  of  product  H2,  as  a 
function  of  the  air  ratio  at  1023  K.  The  cathode  stream  compositions  are  taken 
as  46  mol%  H2/54mol%  H20. 
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Air  ratio 

Fig.  8.  Electrical  energy  consumed  by  the  stack,  per  mole  of  product  H2,  as  a 
function  of  the  air  ratio,  for  average  current  densities  of  7000  and  5000  Am-2, 
and  an  inlet  temperature  of  1023  K. 

ratio,  complements  the  enhanced  convective  cooling  supplied 
by  the  increased  air  flow.  During  endothermic  operation,  on 
the  other  hand,  the  increased  thermal  energy  consumption  com¬ 
petes  against  the  enhanced  convective  heating.  This  is  consistent 
with  the  exothermic  and  endothermic  stack  temperature  distribu¬ 
tions  presented  in  Figs.  5  and  6,  respectively,  in  which  a  greater 
change  in  the  stack  temperature  is  observed  for  the  exothermic 
stack  than  for  the  endothermic  stack  between  the  air  ratio  of  0.4 
and  7. 

Fig.  8  presents  the  electrical  energy  consumed  by  exothermic 
and  endothermic  stack  operation  at  average  current  densities  of 
7000  and  5000  A  m-2,  respectively,  as  a  function  of  the  air  ratio. 
From  the  previous  discussion,  the  electrical  energy  consumed 
by  the  reaction  at  a  given  temperature  is  expected  to  decrease 
with  an  increasing  air  ratio.  Flowever,  the  electrical  energy  con¬ 
sumed  across  the  stack,  here  operating  in  exothermic  mode,  is 
predicted  to  increase.  This  is  due  to  the  increase  in  the  irre¬ 
versible  losses,  associated  with  the  decreased  stack  temperature, 
which  prevails  over  the  decrease  in  the  reversible  potential.  For 
the  endothermic  stack,  on  the  other  hand,  both  the  reversible 
potential  and  irreversible  losses  decrease  with  an  increasing  air 
ratio.  The  overall  reduction  in  the  electrical  energy  consumption 
of  the  endothermic  stack,  shown  in  the  figure,  indicates  that  the 
electrical  efficiency  of  the  stack  is  improved  at  higher  air  ratios. 
During  the  selection  of  the  operating  conditions,  however,  such 
advantages  accruing  from  the  increased  air  ratio  must  be  off¬ 
set  against  the  significant  extra  energy  required  to  generate  the 
increased  air  flow  rate. 

As  previously  mentioned,  SOEC  stacks  are  not  expected  to 
maintain  an  effective  temperature  control  when  operated  near 
thermoneutral  conditions.  The  cathode  and  anode  stream  inlet 
temperature  for  which  the  stack  operates  in  the  thermoneutral 
mode  is  presented  in  Fig.  9,  as  a  function  of  the  average  cur¬ 
rent  density.  The  figure  is  produced  by  introducing  to  the  model 
an  additional  constraint  by  which  the  sum  of  all  the  enthalpy 
changes  involved  in  the  electrolysis  reaction  along  the  stack  is 


Average  current  density  (A  m'2) 

Fig.  9.  Cathode  and  anode  stream  inlet  temperature  corresponding  to  the  ther¬ 
moneutral  stack  operation  as  a  function  of  the  average  current  density.  The  air 
ratio  is  fixed  at  7. 

matched  by  the  electrical  energy  consumed  by  the  stack.  That 
is  to  say  that  the  supplied  electrical  power  not  only  provides  the 
electrical  energy  necessary  for  the  electrolysis  but  also  allows  the 
generation,  through  the  irreversible  losses,  of  the  precise  quan¬ 
tity  of  the  thermal  energy  required  for  the  reaction.  Note  that  the 
curve  in  Fig.  9  represents  the  stack  operation  where  thermoneu¬ 
tral  conditions  are  met  at  every  location  along  the  stack  and 
not  the  situation  in  which  the  thermoneutral  conditions  are  only 
satisfied  at  the  inlet  by  the  local  values  of  the  current  density, 
temperature  and  stream  compositions.  Furthermore,  as  both  the 
thermal  energy  consumed  by  the  electrolysis  reaction  and  the 
heat  produced  by  the  irreversible  losses  depend  on  the  stream 
compositions,  the  thermoneutral  conditions  predicted  here  for 
the  IT-SOEC  stack  with  air-fed  anode  channels  would  be  differ¬ 
ent  to  those  expected  for  the  IT-SOEC  stack  with  pure  O2  anode 
streams. 

The  regions  below  and  above  the  curve  in  Fig.  9  cor¬ 
respond  to  the  conditions  for  exothermic  and  endothermic 
operations,  respectively.  For  average  current  densities  of  7000 
and  5000 Am-2,  an  inlet  temperature  of  1023 K  results  in 
the  stack  operating  near  the  thermoneutral  point.  However, 
by  adjusting  the  inlet  temperature  for  the  7000 Am-2  case 
to  923 K  and  that  for  the  5000 Am-2  case  to  1123K,  the 
stack  can  be  driven  at  more  exothermic  or  endothermic  condi¬ 
tions,  respectively.  The  cathode  stream  temperature  distributions 
along  such  exothermic  and  endothermic  stacks  are  presented  in 
Figs.  10  and  11,  respectively.  The  shift  away  from  the  ther¬ 
moneutral  point  is  evident  in  both  figures  with  an  increased 
temperature  difference  between  the  inlet  and  outlet.  For  an  air 
ratio  of  7,  such  a  temperature  difference  is  as  large  as  88  K  during 
exothermic  operation  and  63  K  during  endothermic  operation. 
The  consequent  increase  in  the  influence  of  the  air  ratio  on  the 
stack  temperature  is  also  observed,  with  the  outlet  temperature 
ranging  by  more  than  120  and  90  K  in  Figs.  10  and  11,  respec¬ 
tively.  The  manipulation  of  the  air  ratio  is  considered  to  provide 
an  effective  means  of  temperature  control  under  these  condi- 
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Dimensionless  axial  position 

Fig.  10.  Cathode  stream  temperature  along  the  stack  for  an  average  current 
density  of  7000  A  m-2,  an  inlet  temperature  of  923  K  and  air  ratios  of  0.4,  7  and 
14. 


7000  Am  2,  an  inlet  temperature  of  1023  K  and  air  ratios  of  0.4,  7  and  14. 


tions.  The  dynamic  response  of  the  exothermic  and  endothermic 
stacks  operated  under  these  conditions  will  be  presented  in  a 
subsequent  publication. 

3.3.  Relationship  between  the  stack  temperature  and 
current  density 

Figs.  12  and  13  show  the  local  current  density  distributions 
for  the  IT-SOEC  stack  with  air-fed  anode  channels  operated  at 
average  current  densities  of  7000  and  5000  A  m-2,  respectively. 
The  air  ratio  is  varied  between  0.4  and  14  while  the  inlet  temper¬ 
ature  is  maintained  at  1023  K.  In  both  figures,  the  curve  displays 
little  dependence  on  the  air  ratio  and  decreases  in  value  towards 
the  outlet.  Generally  along  an  SOEC  stack,  low  local  current 


densities  are  observed  where  increased  voltage  requirements 
exist.  Those  are  the  locations  with  higher  H2  and  O2  partial 
pressure,  lower  FEO  partial  pressure,  and  reduced  temperature. 
Therefore,  as  H2  and  O2  are  generated  and  H2O  is  consumed 
along  the  stack,  the  resulting  change  in  the  stream  composi¬ 
tions  tends  to  decrease  the  local  current  density  towards  the 
outlet.  Furthermore,  the  positive  temperature  gradient  along  an 
exothermic  stack,  and  the  negative  temperature  gradient  along 
an  endothermic  stack,  respectively  act  to  increase  and  decrease 
the  local  current  density  towards  the  outlet.  During  exothermic 
operation,  the  effects  of  the  changing  stream  compositions  and 
the  positive  temperature  gradient  compete  against  one  another. 
Here,  the  decrease  in  the  local  current  density  along  the  stack 
in  Fig.  12  indicates  that  the  effect  of  the  changing  stream  corn- 


density  of  5000  Am  2,  an  inlet  temperature  of  1123  K  and  air  ratios  of  0.4,  7 
and  14. 


Fig.  13.  Local  current  density  along  the  stack  for  an  average  current  density  of 
5000  A  m-2,  an  inlet  temperature  of  1023  K  and  air  ratios  of  0.4,  7  and  14. 
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7000  Am  2,  an  inlet  temperature  of  923  K  and  air  ratios  of  0.4,  7  and  14. 

positions  prevails.  This  is  consistent  with  the  small  temperature 
gradient  observed  in  Fig.  5,  which  is  not  expected  to  significantly 
influence  the  local  current  density  distributions.  For  endother¬ 
mic  operation  in  Fig.  13,  the  changes  in  stream  compositions  and 
the  negative  temperature  gradient  both  contribute  to  decreasing 
the  local  current  density  along  the  stack. 

The  local  current  density  distribution  for  the  IT-SOEC  stack 
with  air-fed  anode  channels  operated  at  an  average  current  den¬ 
sity  of  7000  Am-2  and  an  inlet  stream  temperature  of  923  K 
is  presented  in  Fig.  14,  and  that  for  the  IT-SOEC  stack  with 
air-fed  anode  channels  operated  at  an  average  current  density 
of  5000  Am-2  and  an  inlet  stream  temperature  of  1123K  is 
illustrated  in  Fig.  15.  In  Fig.  14,  the  increase  in  the  local  current 
density  along  the  exothermic  stack  indicates  that  the  effect  of  the 
changing  stream  compositions  is  overwhelmed  by  the  effect  of 
the  temperature  gradient.  Here,  the  local  current  density  distribu- 


Fig.  15.  Local  current  density  along  the  stack  for  an  average  current  density  of 
5000  A  m-2,  an  inlet  temperature  of  1 123  K  and  air  ratios  of  0.4,  7  and  14. 


tion  is  significantly  influenced  by  the  large  temperature  gradient 
shown  in  Fig.  10.  For  the  endothermic  stack,  the  changing  stream 
compositions  and  negative  temperature  gradient  again  both  con¬ 
tribute  to  a  decrease  in  the  local  current  density  towards  the  outlet 
as  shown  in  Fig.  15.  The  increased  dependence  of  the  local  cur¬ 
rent  density  distribution  on  the  air  ratio  in  Figs.  14  and  15  is  a 
consequence  of  the  enhanced  impact  of  the  air  ratio  on  the  stack 
temperature  gradient,  evident  in  Figs.  10  and  11. 

As  the  current  density  is  directly  related  to  the  reaction  rate 
through  Faraday’s  law,  assuming  a  current  efficiency  of  100%, 
the  positive  local  current  density  gradient  in  Fig.  14  corresponds 
to  an  increase  in  the  H2  generation  rate  along  the  stack.  Dur¬ 
ing  exothermic  operation,  such  an  accelerated  reaction  acts  to 
increase  the  rate  of  the  temperature  rise  towards  the  outlet. 
Although  there  are  other  factors  which  influence  the  temper¬ 
ature  distribution  along  the  stack,  this  increase  in  the  reaction 
rate  contributes  to  creating  the  temperature  distribution  to  con¬ 
cave  upwards  in  Fig.  10.  Such  an  effect  is  more  pronounced 
for  operation  with  a  low  air  ratio  where  the  stack  temperature 
is  not  dominated  by  the  convective  heat  transfer  between  the 
cell  components  and  air  flow.  For  endothermic  operation,  on  the 
other  hand,  the  negative  local  current  density  gradient  in  Fig.  15 
and  the  associated  reduction  in  the  reaction  rate  along  the  stack 
tend  to  decrease  the  rate  of  the  temperature  fall  towards  the  out¬ 
let.  Such  an  effect  is  observed  in  Fig.  1 1  as  the  stack  temperature 
distributions,  which  also  concave  upwards. 

4.  Conclusions 

A  one-dimensional  distributed  dynamic  model  of  a  cathode- 
supported  planar  IT-SOEC  stack  with  the  air  flow  introduced 
through  the  cells  has  been  presented.  The  model  has  been 
employed  to  study  the  steady  state  behaviour  of  such  an  SOEC 
and  the  prospect  for  stack  temperature  control  through  varia¬ 
tion  of  the  air  flow  rate.  Steady  state  performance  analysis  has 
shown  that  the  dilution  of  the  product  O2  by  the  air  on  the 
anode  side  of  the  cells  causes  a  reduction  in  the  reversible  cell 
potential.  Activation  overpotentials  were  predicted  to  provide 
the  largest  contribution  to  the  irreversible  losses  while  concen¬ 
tration  overpotentials  remained  insignificant  throughout,  in  spite 
of  the  electrode- supported  nature  of  the  cells.  Evaluation  of  the 
control  strategy  showed  that  an  increase  in  the  air  flow  rate  pro¬ 
vides  enhanced  cooling  and  heating  during  the  exothermic  and 
endothermic  stack  operation  respectively,  prompting  the  tem¬ 
perature  to  become  more  uniformly  distributed  along  the  stack. 
However,  only  a  small  dependence  of  the  temperature  on  the  air 
flow  rate  was  observed  for  a  stack  driven  at  conditions  near 
thermoneutral  operation,  indicating  that  this  operating  mode 
should  be  avoided  from  a  control  perspective.  To  further  analyse 
the  control  strategy,  the  dynamic  response  of  the  stack  will  be 
presented  in  the  future. 
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